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We study the dynamical evolution of cosmological models with the Robertson- Walker symmetry with 
a scalar field non-minimally coupled to gravity and barotropic matter. For this aim we use dynam- 
ical system methods. We have found a type of evolutional path which links between all important 
events during the evolution, the cosmological singularity of finite time, inflation, radiation and matter 
dominating epoch and the accelerated phase expansion of the universe. We point out importance of 
finding the new generic solution called a twister solution for a deeper description of the evolution of the 
Universe. We demonstrate that including the non-minimal coupling leads to a new, richer evolutional 
cosmological scenario in comparison to the case of minimal coupling. 

1 Introduction 

The standard method of description matter content in cosmology bases on the concept of perfect fluid 
approximation, where pressure and energy density satisfy the equation of state. Different cosmological 
epochs constitute solutions of the Einstein equations corresponding different forms of the equation of 
state, usually postulated in a linear form with respect to the energy density. On the other hand if we 
study very early stages evolution of the Universe then characterization of matter content in terms of 
the barotropic equation of state is not adequate. In the quantum epoch including matter in the form of 
scalar field with the potential seems to be more suitable. We propose to describe matter in the form of 
both barotropic matter and single scalar field with the potential. We also assumed that there is present 
a nonzero coupling constant scalar field to the gravity. 

A non-minimal coupling appeared naturally in quantum theory of the scalar field as generated by 
quantum corrections or required by renormalization of the theory^ The value of this coupling constant 
becomes important for cosmology. Recently, it has been constructed an extended model of inflation with 
a non-minimal coupling between the inflaton field and the Ricci scalar curvatureP^HEl was shown in 
particular that the non-minimal inflation can be realized within the Standard Model (SM) or minimal 
extension of itP In the SM the Higgs boson inflation and dark matter is considered by Clark et alP 
The non-m inimal coupling plays also an important role in the context of description of quintessence 
epoch!™ 

The main aim of this paper is to present a generic solution of the non-minimal coupling cosmology 
called the twister solution because of the shape of its trajectory in the phase space. We explore here 
methods of dynamical systems because their advantage of representing all solutions for all admissible 
initial conditions. The phase space is organized by critical points which represent stationary solutions 
for which right-hand sides of the dynamical system vanish and trajectories joining them represent the 
evolution of the system. A new type of evolution appears only if the coupling constant is different 
from minimal and conformal. This solution leads naturally to the quintessence epoch through the twister 



solution. We characterize properties of this acceleration by calculation of so called state- finder parameters. 

We assume the spatially fiat FRW universe filled with the non-minimally coupled scalar field and 
barotropic fluid with the equation of the state coefficient w m . The action is 



S = i / '1 l .'- v -// 
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where k 2 = 8irG, e = +1, —1 corresponds to canonical and phantom scalar fields, respectively, the metric 
signature is (—,+,+,+), R = 6 ^| + is the Ricci scalar, a is the scale factor and a dot denotes 

differentiation with respect to the cosmological time t and U(<f>) is the scalar field potential function. S m 
is the action for the barotropic matter part. 

The dynamical equation for the scalar field we can obtain from the variation SS/Sip = 



<f> + 3H(j> + £R(t> + eU'((f>) = 0, 
and energy conservation condition from the variation SS/Sg^" = 

£ = £ tU 2 + e3£# V + e3£H(<f> 2 y + U(t/>) + p m - ^H 2 . 
Then the Einstein equation for the flat FRW model and the conservation condition read 
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where the energy density and the pressure of the scalar field are 

P4, = e\tf + U{cj))+eUH 2 (j) 2 +e3ZH(<j) 2 l 
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= e~(l - A04> 2 - U{4>) + e^H{4> 2 )- e2£(l - Qi)H4> 2 - e3£(l - 8£)#V + 2^U'{cp). (6) 



2 Non-minimal scalar field cosmology as a dynamical system 



In what follows we introduce the energy phase space variables x = -^jjj, V = > z = ~!jf$i which 

are suggested by the conservation condition 



~\~ q rj2 P ra ^0 1- 



ZH 2 ^ 3H 2 
The acceleration equation can be rewritten to the form 

H 3 
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where the effective equation of the state parameter reads 
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where A - ^ \ dc/ W 
wneie a K uw d<f> . 

The dynamics of the model can be written down as the 4-dimensional autonomous dynamical system 



in variables x, y, z and A where differentiation is with respect to time r defined as -4- 
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is assumed to be the function of 
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A then function z — z{\) can be determined from z(X) = — J[A 2 (r(A) — 1)] _1 dA and the system can be 



further reduced to the form of the 3-dimensional dynamical system 
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if we postulate the exact form of T, for example T(A) = 1 — (a + /3A + 7A 2 ) . Note that this discussion 
is not restricted to the specific potential function but is generic in the sense that it is valid for any 
function T(A) for which z(X) exists. In general there are various potential functions commonly used in 
the literature of the subject^ 

The qualitative analysis of differential equations allows to study all properties of the solutions without 
solving the dynamical system. First, we calculate the critical points which correspond mathematically 
vanishing right-hand sides of the system and physically stationary states. Then we linearize the system 
around these points. The information about the character of critical point and their stability is contained 
in eigenvalues of the linearization matrix. If all real parts of eigenvalues are negative (positive) then 
critical point is stable (unstable). If the all eigenvalues are real of different signs then we have critical 
point of saddle type organized through the stable and unstable submanifolds. It is interesting that for our 
system all important events during the cosmic evolution are represented by these critical points which 
traces a gen eric trajectory representing the evolution of the universe. Let us enumerate all these points 
(see alsoL^. 

1) Finite scale factor singularity; eigenvalues are l\ = 60 I2 = 60 ^3 = 120 an( i critical point are an 
unstable node for positive coupling £ > and a stable node for negative coupling £ < 0. 

2a) Fast-roll inflation; l\ = 0, I2 = 120 ^3 = — 120 an d the critical point is non-hyperbolic. 
2b) Slow-roll inflation; l\ = I2 — I3 — and the critical point is degenerated. 

3) Radiation domination epoch generated by the non-minimal coupling; there are two critical points. 
For the phantom scalar field and £ > eigenvalues are l\ — 0, I2 > 0, I3 < 0. For the canonical scalar 
field and £ > eigenvalues are h — 6£(1 — 3w m ), h — 120 ^3 = — 60 

4) Matter domination epoch; l xft = -| [(1 - w m ) ± [(1 - w m f - ^£(1 - 3W™)] 1 / 2 ] , l 3 = |(1 + w m ), 
and the critical points are non-degenerated for w m 7^ —1 and w m 7^ \. 

5) The present accelerated expansion epoch. 

In the most general case without assuming any specific form of the potential function we are unable 
to find coordinates of this point. In spite of this we are able to formulate general conditions for stability 
of this critical point. From the Routh-Hurwitz criterion we have that the following conditions should be 
fulfilled to assure stability of this critical point 



Re[h,2, 3 ] < ^ 3£^||(y 5 *) 2 > where h(X) = Xz(X) 2 + 4z(X) - 
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To emphasize the acceleration epoch in the twister solution, the state-finder parameters are calculated: 
q = ~d/aH 2 is the deceleration parameter and r = a /aH 3 (see Fig. [IJ. 



3 Conclusions 



In this paper we pointed out the presence of the new interesting solution for the non-minimally coupled 
scalar field cosmology which we called the twister solution (because of the shape of the corresponding- 
trajectory in the phase space). This type of the solution is very interesting because in the phase space 
it represents the 3-dimensional trajectory which interpolates different stages of evolution of the universe, 
namely, the radiation dominated, dust filled and accelerating universe. We are able to find linearized 
solutions around all these intermediate phases, and hence, parameterizations for u> e ff (a) in different epochs 
of the universe history. It is interesting that the presented structure of the phase space is allowed only 



Figure 1: The state-finder parameter diagnostic for the twister quintessence scenario and 2 (A) = — . The trajectories (black 
£ = 6, a = —6, blue £ = 4, a = —4 and purple f = 2, a = —2) represent a twister type solution which interpolates 
between R — the radiation dominated universe (a saddle type critical point), the matter dominated universe (an unstable 
focus critical point) and Q - the accelerating universe (a stable critical point). 



for non-zero value of the coupling constant, therefore it is a specific feature of the non-minimally coupled 
scalar field cosmology. Our general conclusion is that in the description of dynamical complexity the 
cosmic evolution including both barotropic matter and non-minimally scalar field leads to a new richer 
dynamics. 
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